Oscillons in Scalar Field Theories: Applications in Higher Dimensions and Inflation 



O 
O 

(N 



PL, 



> 

oo 



o 



- 1—1 

X 



Marcelo Gleiser 
Department of Physics and Astronomy 
Dartmouth College, Hanover, NH 03755, USA 1 

Abstract 



The basic properties of oscillons - localized, long-lived, time-dependent 
scalar field configurations - are briefly reviewed, including recent results 
demonstrating how their existence depends on the dimensionality of space- 
time. Their role on the dynamics of phase transitions is discussed, and 
■ it is shown that oscillons may greatly accelerate the decay of metastable 

' vacuum states. This mechanism for vacuum decay - resonant nucleation 

- is then applied to cosmological inflation. A new inflationary model is 



1 proposed which terminates with fast bubble nucleation. 



1 Introduction 



During the past three decades, static, spatially-localized field configurations 
have been of great interest in relativistic field theories ^ |2] • These nonpertur- 
£SJ . bative solutions of nonlinear field equations have been shown to exist in a wide 

variety of models, with or without nontrivial vacuum topology. Configurations 
with stability guaranteed by vacuum topology are called topological solitons. 
Examples include kinks, strings, monopoles, and textures [3]. When the config- 
uration's stability comes from conserved charges, they are called Q-balls 0] or 
nontopological solitons Both topological and nontopological solitons have 
many applications in high-energy particle physics and cosmology Q][3J. Their 
widespread appeal relies on two key properties: first, they are stable, that is, 
they maintain their spatial profile through time. (Of course, if the configurations 
are allowed to scatter they may or not remain stable and form bound states.) 
Second, their energy density is spatially localized. Thus, they have been used 
$— i ' to model hadrons QJE], as possible signatures of primordial phase transitions, 

or to generate the large-scale structure of the Universe PJ. 

The purpose of the present work is to present results pertaining to a related 
but, in this author's view, equally important class of spatially-localized field 
configurations. Their key difference from the solitons mentioned above con- 
cerns their temporal behavior: they are time- dependent, as opposed to static, 
solutions of the equations of motion. What makes them potentially interesting 
for many applications is that, even though they will eventually radiate their 
energy to spatial infinity, the process is quite slow. If their lifetime is longer 
than the typical time-scale in the system they will, for all practical purposes, 
behave as solitons. The remarkable fact about these long-lived, time-dependent 
configurations - called oscillons - is that they owe their longevity to the nonlin- 
earities in the system: no globally-conserved charges or topologically nontrivial 
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boundary conditions are needed (although they may help extend the oscillon's 
lifetime). As such, they can be found in a much broader class of models. 

In the next section oscillons are briefly introduced, followed by an exploration 
of their properties in an arbitrary number of spatial dimensions 5' ■ In section 
3, the possible effect of oscillons on first order phase transitions is discussed 
UJ. Finally, in section 4, they are applied to cosmology, in particular to a new 
two-field model of inflation that ends with rapid bubble nucleation. 

2 Oscillons in d Dimensions: Basics 

Oscillons were first shown to exist in the context of simple scalar field theo- 
ries with symmetric and asymmetric double- well potentials |S] . They were 
found to be spherically-symmetric, time-dependent solutions of the equations of 
motion obtained from the c£-dimensional LagrangianjS], 

where a dot means time derivative. The c?-dimensional spatial volume element 
can be written as d d x = Cdr( d ~ 1 * , dr, where ca = 2it d / 2 /T{d/2) is the surface 
area of a <i-dimensional sphere of unit radius. An oscillon is a time-dependent, 
spatially-localized solution of the equation of motion <f> — \/ 2 d 4> — —V, character- 
ized by a nearly constant energy, Ease , and by a persistent nonlinear oscillation 
about its core at r = 0. Think of a rubber sheet that is pinched and then let 
go. An oscillon would be a localized oscillation on the sheet that doesn't get 
radiated away for a long time. To find an oscillon, write the scalar field as 

0(4, r) = [cj> c (t) - <j> v ] exp[-r 2 /R 2 } + cj> v , (2) 

where 4> c (t) is the core value of the field [<fi(t,r = 0)], and <j) v is its asymptotic 
value at spatial infinity, determined by V(<p). Using the Gaussian profile of eq. 
[21 it has been shown that oscillons are present whenever the energy of the initial 
configuration (j)(0, r) is larger than E osc and R > R OS c, where R osc is the oscillon 
radius that can be obtained analytically 7 . In addition, the potential must 
have a portion where V" < and the field must probe this portion at r ~ 0. In 
fact, the oscillon owes its stability to this "spinodal instability" in the potential 

m 

Writing A(t) — (f> c [t) — 4> v and parameterizing the potential as 5 

h 

3=1 J ' 

where the g^s are constants and the vacuum energy V{<f> v ) is subtracted from 
the potential to avoid spurious divergences upon spatial integration, one obtains 
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the equation of motion 
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Expanding the amplitude as A(t) = + SA(t) and linearizing the eom 

with 5 A ~ e lwt , it is found that oscillons may only exist if w 2 < 0. For quartic 
potentials with > (that is, symmetric and asymmetric double wells), this 
conditions implies 

R 2 > ~r A ~ ■ (5) 



i( 



2 3 / 2 
3 



(V'") 2 
yiv 



v 



That is, oscillons only exist for lumps above a certain critical size. For d = 3, 
i? osc ~ 2.42, while for d = 6, R osc ~ 7.5. Note that since the denominator of 
cq. [S] must be positive, there is also an upper critical dimension for oscillons to 
exist, 



d < Int 



] n i V"V 1V . 



In 
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For a symmetric double- well, d c = 6 p]. 

The fact that for these simple models oscillons can only exist below a certain 
number of spatial dimensions raises an interesting possibility. Let's assume that 
this model carries the information of the scalar sector of a more realistic field 
theory. (Far hi et al. recently found oscillons in a SU(2) model 0.] Then, if 
extra dimensions exist and are large it is conceivable that oscillons could 
be produced in collisions with high enough energy, E co \ > E osc . Since their 
lifetime is much longer than the typical time-scales in the system, they would 
appear as a late fireball concentrated in a volume ~ Ro SC - Furthermore, since 
E sc ~ ^(n/2) d / 2 d d ~ 1 , finding an oscillon would provide information about the 
dimensionality of spacetime [S]- In other words, oscillons could be used as probes 
for the existence of higher dimensions. 



3 Resonant Nucleation 

Recently, oscillons have been shown to emerge dynamically when a system is 
quenched from a single to a double well potential This means that, under 
certain conditions, they are spontaneously produced during the noncquilibrium 
dynamics of the system. [The same results would be obtained at T = 0, as 
long as the field starts in, say, a Gaussian superposition of momentum modes.] 
The system was initially prepared in a thermal state on a single well and then 
quenched to a symmetric double well, as is common in Ginzbug-Landau models 
of phase transitions ^2] . The key point is that the system was prepared so that 
after the quench the field was still localized in one of the wells. For example, the 
potential could have been changed from V(<p) = ((/> + l) 2 to V(<p) — \{4> 2 — l) 2 - 
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Figure 1: The number of oscillons nucleated between t and t + St at T = 0.2, 
with St = 1. The global emergence is evident early in the simulations. Inset: 
the probability distribution of radii and periods of oscillation for the oscillons 
nucleated. 



The field would then stay trapped around the <fi = — 1 minimum. This is what 
happens, for example, during the initial stages of a first order phase transition, 
when the field is initially in a metastable state. More generally, one can write a 
potential 

V{4>) = 2> - ^ + ^ , (7) 

where a is a tunable interaction. Initially, a — and the field is in a single well. 
Then a acquires a positive nonzero value and the potential becomes a double 
well. For a — 3/2 the double well is symmetric; for a > 3/2 the minimum at 
(j) = is metastable jllj . 

The net effect of the quenching is to decrease the effective mass of the field 
around the minimum. This key point may not be apparent at tree level, es- 
pecially around = 0. However, one must remember that the mass must be 
corrected due to thermal and/or quantum fluctuations. One way of incorpo- 
rating these changes is to adopt a simple Hartree approximation, shifting the 
mass from m 2 to m 2 H — m 2 + |(^ 2 ). Before the quench, the Hartree potential 
is (upon rescaling the couplings) Vh{4>) = |[1 + |(0 2 )]]</ )2 - After the quench, 

1 a 1 

V H (<P) = [1 - m 2 H ] cj> + -m 2 H cf> 2 - -0 3 + -0 4 . (8) 

Thus, the quench induces a shift in the minimum from <p — to ~ |(0 2 ). 
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Figure 2: The evolution of the order parameter (f> ave (t) at T = 0.22 for 
several values of the asymmetry in d = 2. (From left to right, a = 
1.746,1.56,1.542,1.53,1.524,1.521,1.518. Each curve is an ensemble average 
over 100 runs. The inset shows the effective Hartree potential for the same 
values of a. 

This shift, plus the fact that Vh(4>) acquires a negative derivative at the origin, 
explains what happens: the field's zero mode starts oscillating about the shifted 
minimum, with an amplitude controlled both by a and T. These oscillations 
may induce parametric amplification of a band of A:-modes which, for a range 
of parameters, trigger the emergence of oscillons 11 . The oscillations trans- 
fer their energy to higher momentum modes via nonlinear scattering and the 
oscillons eventually disappear as the system reaches equipartition. 

Fig. 1 shows the number of oscillons nucleated as a function of time. Two 
points are important: first, oscillons emerge in synchrony. Second, this initial 
synchrony is gradually lost as equipartition is achieved. 

What happens if the potential is asymmetric? In this case, one is concerned 
with the way a quenched field will decay from a metastable (or false vacuum) 
state to the global minimum. As is well known, if the field is initially well- 
localized about the false vacuum state, it will decay by nucleating bubbles of a 
critical size with nucleation rate per unit volume T ~ T 4 exp[— E(T)/T], where 
T is the temperature and E(T) is the energy of the critical bubble or bounce 
configuration [E]. At T = 0, the quantum decay rate is Y ~ M 4 exp[— Sb], 
where M is the typical mass scale and Sb is the d + 1-dimensional Euclidean 
bounce action. 

As in the case of a symmetric double well, the fast quench will generate a 
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distribution of oscillons for a certain range of temperatures. The efficiency of 
this production mechanism is sensitive to the quenching time scale, r quonc h. If 
^quench 3> To, where tq is the equilibration time scale of the field's zero- mode, 
oscillons will be produced. In Fig. 2 we show the variation of the average value 
of 4> after the quench. It can be seen that close to degeneracy (a ~ 3/2) the field 
spends a long time oscillating about the metastable minimum before decaying 
to the ground state. 

The key point is that a fast enough quench will greatly affect the decay rate 
UJ. As shown in Ref. 0, a fast quench will promote a much faster decay: for 
a range of parameters, the decay rate goes from an exponential to a power law 
suppression, 

r ~ T 4 [E{T)/T]- B , (9) 

where B can be obtained numerically. For d = 2, it was shown that B ~ 
2.9 ± 0.46 for the range of temperatures where oscillons are copiously produced 
6 . In Fig. 3 the decay time-scale is plotted against the bounce energy for 
different temperatures. It is clear that a power law is an excellent approximation. 

Work in d = 3 is currently under way. Preliminary results indicate that the 
power law behavior remains, with B ~ 2. 

The main consequence of the above result is that whenever the quench- 
ing happens fast, using an exponentially-suppressed decay rate as dictated by 
false-vacuum decay theory is simply incorrect. Since T ~ exp[—E(T)/T] (or 
~ exp[— Sy at T = 0) is widely used in applications ranging from early universe 
cosmology (in particular at GUT scales, where H is fast) to models of quark- 
hadron phase transitions in collisions such as in RHIC, one should be careful 
about their range of validity. It may very well be that transitions that were 
considered slow in the past are actually quite fast. As an application of what 
was called "resonant nucleation" (RN) in Ref. UJ, I will briefly investigate its 
possible effects on inflation. 

4 Resonant Inflation: Can Old Inflation be Res- 
cued? 

The simple elegance of the original "Old" Inflation (01) scenario proposed by 
Guth in 1981 has, since then, inspired many variations More than just 

the elegance of its formulation, based on a single scalar field decaying from an 
initial metastable state to a lower-energy state by bubble nucleation, the original 
01 model had a clear connection with particle physics: the inflaton was to be the 
same scalar field promoting the symmetry breaking of Grand Unified models, 
linking early-Universe cosmology to high-energy particle physics. In fact, it is 
this particle physics connection that motivated and motivates the widespread 
use of scalar fields in early-Universe physics. 

Unfortunately, Guth's original proposal didn't work. As he himself argued, 
and then Linde, and Albrecht and Steinhardt ^B], the bubble- nucleation rate 
could not compete with the exponential expansion rate of the Universe: the 
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Figure 3: The decay time-scale trn for resonant nucleation as a function of 
critical nucleation free-energy barrier E(T)/T at T = 0.18,0.2, and T = 0.22. 
The best fit is a power-law with exponent B ~ 2.44 for T = 0.22, and B ~ 3.36 
for T = 0.18. 

transition would never end. Roughly, while bubble walls expanded with the 
speed of light, their centers receded from each other exponentially fast, making 
it impossible for the walls to touch, the bubbles to coalesce, and the transition 
to complete. Old Inflation gave rise to a universe with inhomogeneities incom- 
patible with the observed smoothness of the cosmic microwave background |17| . 
Guth and Weinberg [TQ, an( i later Turner, Weinberg, and Widrow [TJJ, per- 
formed a detailed analysis of the constraints needed to render OI and OI- inspired 
scenarios viable. They concluded that a strong (or, equivalently, slow) first or- 
der phase transition based on a single scalar field could not be made to work: 
the ratio of decay rate to the expansion rate per unit volume [H A ~ (T 2 /Mpi) 4 ], 
had to be sufficiently small 



initially, so that early bubbles didn't produce inhomogeneities during nucleosyn- 
thesis and on the CMB. (For T ~ T A exp[-E(T)/T] and T GUT = 10 15 GeV, this 
implies that E(T)/T > 46.1 initially). On the other hand, it had also to grow by 
the end of inflation (e — > 97r/4) to guarantee that the transition was completed 
|19|. [This implies E(T)/T < 34.9.] In other words, successful inflation forced 
the decay rate to be time-dependent: small at the beginning of inflation and of 
order unity at the end. As further work has shown, this could be achieved by 
invoking more fields |20| and/or a nonminimal gravitational coupling 21 . 
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Figure 4: Comparison of homogeneous nucleation (HN) and resonant nucleation 
(RN) with power B in an expanding Universe. For a fixed B and nucleation 
barrier (3 — E(T)/T (or Sb at T = 0), the line denotes equality. Values above 
the curve imply faster HN, while those below imply faster RN. 

Given what we have learned in the previous section about resonant nucle- 
ation, it is natural to wonder whether such effects can play a role on inflation. 
If we write ern — T 4 exp[— E(T)/T] to represent the ratio of eq. ^3 using the 
homogeneous nucleation rate, and ern — T 4 {E(T)/T]~ B the ratio using the 
RN rate, equality is attained whenever 



where (3 = E{T)/T (or = Sb at T = 0). In Fig. 4 B is shown for representative 
values of the nucleation barrier j3. The line denotes Ehn/ern = 1- The squares 
denote the limits imposed by the inflationary constraints of Ref. JE] . For these 
values of /3, unless B > 9, which is very unlikely, resonant nucleation rates are 
always faster. In d — 2, where B ~ 3, it is clear that ehn/srn < 1 for all 
realistic values of /?, not a surprising result. 

Why is this useful for inflation? For successful inflation with HN, the con- 
straints of Ref. limit the nucleation barrier f3 to be fairly small 0(~ 40). 
[See Fig. 21] However, calculations of bounce actions show that f3 usually scales 
with inverse powers of coupling constants. These two requirements compete 
with each other, making it hard to have small nucleation barriers with small 
couplings. Applying the percolation constraint to the RN rate, one obtains 
(3 s ~ 10 16 . For B = 3, this gives /3 ~ 10 16 / 3 : small couplings (or, equivalently, 
large barriers) are easier to accommodate with RN, the first reason why it may 
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be useful for inflation. 

The second and most important reason is that RN makes it much easier 
to complete the transition. Clearly, if some mechanism capable of producing 
the same net effect as the fast quenching responsible for RN was present in 
the early Universe, an initially slow first order transition could become fast 
at some point, going from an exponential to a power law decay. In this way, 
even a potential with a large initial barrier would not be an impediment to the 
successful termination of inflation. One possible way of implementing RN in 
cosmology is to invoke a second field ip that couples to the nucleating field (p in 
a way somewhat reminiscent of hybrid inflation. In that model the inflaton <f> 
is coupled quadratically to another scalar field ip which has a symmetric double 
well potential [2"U|: 

V(</>, fl>) = ^ (M 2 - A-0 2 ) 2 + ^> + |W • (12) 

Inflation is driven by the energy density V((p, 0) while the inflaton (cp) is rolling 
down along the ip = valley (T£|- As <p reaches a critical value, ip becomes 
spinodally unstable and quickly rolls to one of the minima (or both, but this 
creates domain walls, another problem), terminating inflation abruptly. 

The key difference with the mechanism being proposed here is that bubble 
nucleation still occurs at the end of inflation. A possible name is hence resonant 
inflation (RI): it blends 01 with the physics of resonant nucleation. 

Modify the potential for the field <p that gives rise to RN (e.g. eq. 0) by 
coupling another field (ip) quadratically to it as follows, 



V(4>, 1>) = \ (m 2 + 5 V) - |0 3 + ^ 4 + \m% ^ + \V(ct> + ,0)\ , (13) 



where <p+ is the value of <p at the global minimum of V(cp, 0) so that V(0, ip) 
provides the net vacuum energy responsible for inflation. [Note that here the 
inflaton is ip.) Inflation lasts while ip is rolling down the </> = valley. Notice that 
the mass term for <p, AfJ = m 2 + g 2 ip 2 , decreases as ip rolls down its potential. 
While M 2 > a 2 /2A, the only minimum in the cp direction is at (p — 0. However, 
as ip decreases, Af? will eventually drop below a 2 /2A and a new minimum will 



appear at 



At Af? = the two minima 



= f [l + (l-2A/ 2 A/a 2 ) 1 /2 

are degenerate. At this point, as ip continues to approach zero, the minimum 
at (p — becomes metastable. Oscillations in cp, induced by the decrease in 
its mass, will induce RN. This will be true as long as the decrease in Af^, 
M$ ~ B-ipifi, is fast enough. [It was assumed for simplicity that g 2 ip 2 /m 2 <^ 1 
which is not true for very small ip.] 

For RI to work, M^/M^ < H during inflation and M^/M^ > H after it. 

During inflation, with a slow-roll approximation, ipip ~ —-^-ip 2 . We then ob- 



tain, 
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Also, if N is the number of e-folds, ip 2 — ipf g^ 1 -^' wnere V'i(e) is the value 

of the field ip at the beginning (end) of the inflationary period. [For simplicity, 
it was assumed that during inflation ^m^ip 2 > \V(<p+,0)\, that is, inflation is 
dominated initially by the vacuum energy of the inflaton field ip.] Slow-roll 
ends when ip 2 < Mp,/12%. Using this result and eq. El the slow variation 
of M implies, £? 2 (12tt) 1 / 2 < {m/M P {f. [For m ~ 10 16 GeV, g < 4 x 1CT 4 .] 
This condition is also consistent with the approximation g 2 ip 2 /m 2 <C 1 for 
ipi > ip > ip e , that is, during inflation. 

If slow-roll ends when the minimum in <p + appears, we obtain (this is similar 
to the critical condition in hybrid inflation 20 ), 

where we defined for convenience a = m«o. The condition for slow variation 
of during inflation forces the second term on the rhs of eq. Elt° be very 
small. Thus, if we want to impose that the <p+ minimum appears close to the 
end of inflation, we must have «q/2A ~ 1, not a difficult condition to satisfy. 

As inflation ends, ip will start rolling down fast towards the ip = minimum 
and oscillate around it. Since in this regime, M^/M^, ~ {g 2 / m 2 )tpijj , the rapid 
motion of ip will induce the time-dependence in needed to trigger resonant 
bubble nucleation. In order for the transition to end successfully, the percolation 
constraint £rn > 9/47T, must be satisfied. This implies, 

x b „ 4vr (Mpi^ 4 



If m ~ 10 16 GeV and B ~ 2 (as indicated by preliminary results in d = 3), RI 
terminates efficiently if St, < 10 6 . Since the inflationary phase is due to the slow- 
roll dynamics of the ip field and not by the metastable field <p, the percolation 
constraint can be satisfied by a wide range of couplings. Also, since = only 
becomes metastable after the end of slow roll, there is no need to impose the big 
bubble constraint: once ip starts rolling fast at the end of inflation, RN will ensue 
and rapid bubble nucleation and coalescence will quickly reheat the Universe. 
Although several details remain to be worked out, this preliminary analysis 
indicates that resonant nucleation can be successfully applied to inflationary 
cosmology. 
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